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ABSTRACT 

We report Hubble Space Telescope imaging obtained 155 days after the 2006 
outburst of RS Ophiuchi. We detect extended emission in both [Om] A5007A 
and [Ne v] A3426 A lines. In both lines, the remnant has a double ring structure. 
The E-W orientation and total extent of these structures (580 ± 50 AU at d — 
1.6kpc) is consistent with that expected due to expansion of emitting regions 
imaged earlier in the outburst at radio wavelengths. Expansion at high velocity 
appears to have been roughly constant in the E-W direction (v exp = 3200 ± 300 
km s -1 in the plane of the sky), with tentative evidence of deceleration N-S. 
We present a bipolar model of the remnant whose inclination is consistent with 
that of the central binary. The true expansion velocities of the polar components 
are then v = 5600 ±1100 km s _1 . We suggest that the bipolar morphology of 
the remnant results from interaction of the outburst ejecta with a circumstellar 
medium that is significantly denser in the equatorial regions of the binary than at 
the poles. This is also consistent with observations of shock evolution in the X-ray 
and the possible presence of dust in the infrared. Furthermore, it is in line with 
models of the shaping of planetary nebulae with close binary central systems, and 
also with recent observations relating to the progenitors of Type la supernovae, 
for which recurrent novae are a proposed candidate. Our observations also reveal 
more extended structures to the S and E of the remnant whose possible origin is 
briefly discussed. 

Subject headings: binaries: symbiotic - novae, cataclysmic variables - planetary 
nebulae: general - stars: individual (RS Oph) - supernovae: general - white 
dwarfs 



1. Introduction 



RS Ophiuchi is a symbiotic recurrent nova (R N) which had previously undergone re corded 
outbursts in 1898, 1933, 1958, 1967 and 1985 (see Rosinol ll987t iRosino & Iiiimalll987h. with 
possi ble additional outbursts in 1907 (jSchaeferl 12004 ) and 1945 (jOppenheimer fc Mattei 
19931 ). On 2006 February 12.83UT it was observed to be undergoing a further eruption 
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( IHirosawall2006l ) . reaching magnitude V=4.5 at this time. For the purposes of this paper, we 
define this as t = 0. The op tical light cu rve then began a rapid decline, very similar to that 
seen m previous ouib^ eH AAVSC0) - Several lines of evidence are consistent 



with a distance to RS Oph of 1.6 ± 0.3 kpc (lBoddll987l ) 



The RS Oph binary system comprises a red giant star in a 455.72 ± 0.83 day orbit with 



a white dwarf (WD) of mass near t he Chandrasekhar limit (see iDobrzycka fc Kenyonl 11994 



Shore et al. 


1996; 


Fekel et al. 


2000) 



onto the WD surface leads to the conditions for a thermonuclear runaway (TNR) in a similar 
fashion to that for classical novae (CNe). The much shorter inter-outburst period for this 
type of RN compared to CNe is thou ght to be due to a combination of the high WD mass and 



a supposed high accretion rate (e.g. IStarrfield et al.lll985c lYaron et al.ll2005l ). Such models 



lead to the ejection of somewhat lower masses at higher velocities than those for models of 
CN (typically 5000 km s -1 and 10 -8 — 1O~ 6 M respectively for RNe). Spectroscopy of RS 
Oph has indeed shown Ha line emission w ith FWHM = 3930 km s _1 and FWZI = 7540 km 
s- 1 on 2006 February 14.2 {t = 1.37 davs. iBuiJlioO^ . 



Unlike CNe, where the mass donor is a low-mass main-sequence star, the presence of 
the red giant in the RS Oph system means that the high velocity ejecta run into a dense 
circumstellar medium in the form of the red giant wind, setting up a shock system with gas 
temperatures ~ 2.2 x 10 8 K for vg = 4000 km s _1 , where vs is the velocity of the fo rward shock 



runnin g int o the pre-existing wind . X-ray observations of the 2006 outburst by iBode et al. 



( 2006 ) and ISokoloski et all (120061) have confirmed the basic shock model of IBode fc Kahn 
(Il985l ) and lO'Brien et all <jl992h in which RS Oph evolves like a supernova remnant (SNR), 
but on timescales around 10 5 times faster. 



VLBI imaging at 6cm and 18cm began 13.8 days after outburst (jO'Brien et al.ll2006l ). 
The initial image showed a partial ring of non-thermal (synchrotron) emission, of radius 13.8 
AU at d — 1.6 kpc, consistent with emission from the forward shock. Subsequently, more 
extended lobes aligned E-W gradually emerged, with that in the East appearing first. This 
morphology was consistent with that derived f r om m ore rudimentary VLBI observations 77 
days after the 1985 outburst by iTaylor et al.l (119891 ). apparently showing jets of emission. 
O'Brien et al.l (120061 ) proposed a simple geometrical model comprising an expanding double- 
lobed structure, with the major axis perpendicular to the proposed plane of the central 
binary orbit. 

Here we report Hubble Space Telescope (HST) observations of the expanding nebular 
remnant taken at t — 155d and compare them to the structures seen earlier in the radio. We 



1 http:/ /www. aavso.org 
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then go on to explore remnant geometry by comparison to a bipolar model and thus deduce 
more about the circumstances of the outburst and the wider implications of these results. 



2. HST Observations and data reduction 

RS Oph was observed (Prop. ID 11004) in Director's Discretionary Time 155 days after 
outburst on 2006 July 17 using the Advanced Camera for Surveys (ACS) on the HST. The 
0.025"pixel _1 High-Resolution Channel (HRC) was used for the observations combined with 
three narrow-band filters to isolate the Ha+[Nll] (filter F658N), [Om] A5007A (F502N) 
and [Nev] A3426A (F343N) nebular emission lines. In order to reduce readout overheads 
and maximise on-object time, the 512x512 square pixel subarray was used throughout. 

One orbit was spent imaging RS Oph, while a second was used to observe HD 166215, a 
bright nearby standard star of comparable spectral type (V = 9.32, K0). With limited time 
for the RS Oph observations, a balance was struck between deep imaging and dithering to 
increase resolution. Care was taken when planning the exposures not to saturate either the 
primary target or standard star images in any filter. The resulting images were checked and 
verified as unsaturated in every case. Observations of RS Oph were repeated with a 3" offset 
in both RA and declination in order to identify and eliminate ghosts if present. A two point 
dither pattern was used for each observation. Table 1 summarises the resulting observations 
of each target. 

All data were reprocessed using standard procedures outlined in the ACS Data Hand- 
book^! and the Pydrizzk^l and Multidrizzl^l Handbooks, using optimal input parameters to 
maximise the signal to noise ratio of each image. Although extended emission was clearly 
visible in the unprocessed [O ill] images, a number of data analysis techniques were used in 
order to explore in detail the presence of extended emission in the images for each filter. 



PSF profiles generated by TinyTim (IKristl Il995l ) were then subtracted from the RS 



Oph data in each of the emission lines. The PSF was scaled to an extent such that when 
subtracted, it brought the formerly peak pixel value to approximately the same level as those 
surrounding it. After subtraction, the total counts over an area enclosing the central star 
(as a percentage of the pre-subtracted values) are 92% ([Om] and [Nev]) and 65% (Ha). 
It is likely therefore that there is some residual stellar flux in the central regions of the 



2 http: / /www. stsci.edu/hst/acs/documents/handbooks/DataHandbookv4/ ACS Jongdhbcover.html 
3 http: / / stsdas.stsci.edu/pydrizzle 
4 http://stsdas. stsci.edu/pydrizzle/multidrizzlc 
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PSF-subtracted images. 

Using the profile of HD 166215 as the PSF, and also that generated by TinyTim, de- 
convolution using the Lucy- Richardson method was performed on the RS Oph data for each 
emission line. Tests using both the CLEAN and Maximum Entropy techniques produced 
similar results. However, for [Nev] the lower signal-to-noise of the stellar PSF meant that 
only the TinyTim PSF could be used effectively in the deconvolution. 

As part of this study, a re-analysis was carried out of pre-outburst WFPC2 observations 
of the [Om] A5007A line through the 502N filter on 2000 June 12 (Prop. ID 8332). No 
extended emission, even at very faint levels, was detectable. This confirms the results of the 



analysis of these data undertaken by lBrocksopp et al.l (120031 ). 



3. Results and Discussion 

As can be seen in Figured], extended structure was clearly visible in the [O ill] A5007A 
line in the PSF-subtracted image (and was indeed visible in the raw image). Deconvolution 
revealed more detailed structure in both [O ill] and [Ne v] A3426 A. There was also a hint 
of possible extended emission close to the central source in the Ha line, but this was not 
present at a significant level. 

In the deconvolved [O ill] and [Ne v] images, the most striking feature is a double ring 
structure with major axis lying E-W and total (peak-to-peak) extent 360 ±30 mas (580 ±50 
AU at d = 1.6 kp c). The most extended structures seen in the radio (the outer lobes; 



O'Brien et al.l 120061 ) lie along this axis. Assuming ejection at the time of the outburst, the 
expansion rate (from the center) of the optical emission along this axis is 1.2 ± 0.1 mas d _1 
(equivalent to v exp = 3200 ± 300 km s _1 in the plane of the sky). The optical emission is 
also detectable above background in the deconvolved images to a total extent of 520 ± 50 
mas, corresponding to an expansion rate of 1.7 ± 0.2 mas d _1 . We may compare this to 1.4 
±0.3 mas d" 1 for the E-W lobes seen in the radio, taken over 4 epochs from t = 21.5 to 62.7 
days in the 2006 outburst (O'Brien et al., in preparation) and 1.3 mas d -1 for the eq uivalent 



features derived from VLBI observations in the 1985 outburst (ITaylor et al.lll989l ). Thus 
there is evidence that the bipolar emission seen here and in the radio arises from the same 
regions of the remnant, if the expansion velocities in the E-W direction are roughly constant 
over the first 155 days. 

Tracing the ring-like structure of the optical remnant suggests a N-S (peak-to-peak) 
extent of 150±25 mas. This corresponds to an expansion rate of 0.48±0.08 mas d -1 compared 



with 0.62 ±0.01 mas d 1 derived in the radio for the central ring on day 13.8 (jO'Brien et al. 
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20061 ). As noted in lO'Brien et al.l (2006). the form al error on this radio expansion rate is 
likely to be an underestimate; iRupen et al.l (120071 ) also derived radio expansion rates for 
the ring ~ 20% higher. Thus comparison of the HST and radio data tentatively suggests 
deceleration in the N-S direction, but verification awaits further optical observations. 

As noted above, both [Olll] and [Nev] show similar structure. In the case of both 
lines, emission might arise either from regions downstream of the forward shock in the 
red giant wind or from a "precursor" ionised region ahead of the shock, or both (e.g. 
Dopita fc Sutherlandlll996l ). The more marginal presen ce of Ha+[Nn] may be a d irect con- 
sequence of particular shock velocities (again see e.g. iDopita &: Sutherland! Il996l ) coupled 
with decreasing contrast of extended emission with the central star at longer wavelengths. 
It is also of interest to not e that the assump tion from Spitzer observations of infrared fine 
structure and coronal lines (jEvans et al.ll2007l ) of an origin of these lines in the same remnant 
regions is consistent with our observations. 



O'Brien et al.l (120061 ) proposed a simple model for the radio emission comprising a bipo- 



lar structure where the evolution of the optical depth due to free-free absorption led to 
the gradual uncovering of various features, most notably this explained the emergence of 
the outermost radio lobes. Here, we have modeled the optical emission seen in our HST 
images with a "pe a nut-s haped" bipolar structure using the modeling code described in 
Harman &: O'Brien! (120031 ) (see Figure [2]). The shell has finite thickness (0.05 arcsec) and the 
surface brightness distribution is just determined by the line-of-sight path length through 
the nebula. The axial ratio is set at 3:1, consistent with a constant velocity of expansion 
along the major axis and deceleration along the "waist", as implied by the earlier radio 
observations. 

As can be seen from Figure [21 the model reproduces the morphology of the optical 
emission extremely well for an inclination % = 35°, consistent with the major axis lying normal 
to the plane of the binary orbit, where i = 3 0° — 40° was g i yen b y iDobrzycka &: Kenyon 
(119941 ) from radial velocity studies. Following lO'Brien et al.l (120061 ). we also note that the 
2006 outburst occurred at approximately the same binary phase as in 1985, and hence the 
similarity between the rate of expansion of the radio lobes on the plane of the sky in the two 
outbursts would naturally follow. 

We now turn our attention briefly to the features to the E and S of the main remnant 
(Figured]). The eastern "arc" , which is apparent only in the [O ill] image, has no counterpart 
detectable to the W and its closest approach to the centre of the bipolar emission is 0.44±0.01 
arcsec. Similarly, the southern "blob" (also detected in Ha) lies at 1.52 ± 0.02 arcsec from 
the centre. If these features arose in the 2006 outburst, for d — 1.6 kpc, they would therefore 
imply v exp ~ 8000 and ~ 27000 km s _1 respectively, values far higher than those otherwise 
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observed or expected in such an outburst. On the other hand, at this time, one might expect 
light echoes in the surrounding medium to have radial extent far greater than this (~ 16.3 
arcsec, for a scattering medium in the plane of the sky at the distance of the nova). Another 
possibility is that these features are in fact associated with material piling up from previous 
outbursts. Future HST observations will help to clarify their origin, and also explore the 
late-time evolution of the remnant. 



4. Conclusions and future work 



We have detected, in at least two emission lines, the resolved optical counterparts of the 
expanding emission first seen much earlier in the outburst at radio wavelengths. The central 
remnant is then modeled in terms of a bipolar s tructure wh ose inclination is consistent with 
that of the binary orbit. For d — 1.6 ± 0.3 kpc (jBodelll987l ). % = 35°, an expansion velocity, 
v = 5600 ± 1100 km s _1 and v ra d = 4600 ± 900 km s _1 are implied for the material at the 
poles (the formal errors being dominated by the uncertainty in distance here). Unless the 
ejection of material at the explosion site is highly anisotropic (which appears unlikely) this 
reinforces the notion that remnant shaping occurs due to the interaction of the ejecta with 
the pre-existing circumstellar environment. The implication would then be that the densest 
parts of the red giant wind lie in the equatorial regions along the plane of the binary orbit. 
This in turn leads to more marked deceleration of material in this direction compared to the 
poles. The existence of a dense region of the wind, which has not been fully traversed by a 
forward shock (even at these relatively late times) may well explain the apparent presence of 
dust in the circumstellar environment throughout the outburst (Evans et al., in preparation). 
It is from this dense equatorial region that we now suspect most of the early X-ray emission 
also arises as the forward shock driven into this region by t h e out burst ejecta decelerates. 
The rapid evolution of the remnant derived by iBode et al.l (120061 ) is consistent with this 
proposal. 

These observations have wider relevan ce to our unde r stand ing of the shaping of nova 
shells and planetary nebulae For example, iMitchell et al.l (120071 ) have shown that the PN 
Abell 63 comprises narrow lobes which form a "tube" at the ends of which are "ca ps" of emis 
sion . At th e centre of Abe l l 63 is the eclipsing close binary UU Sge (p = 11.2 hrs lBond et al. 



19781 ). and IMitchell et al.l (120071 ) have demonstrated that the inclination of the lobes is that 
of the plane of the central binary. They conclude therefore that the PN was shaped by in- 
teraction during its early evolution with a circumbinary environment with enhanced density 
in the equatorial plane. This is consistent with the results of hydrodynamic modelling of th e 
evolution of the environment of such systems performed by iMastrodemos fe Morris! (119991 ). 
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Although they did not extend their study to semi-detached (rather than detached) systems, 
these authors found that equatorial enhancements in the density of the winds from red giants 
and AGB stars in binary systems were most pronounced for those with the smallest binary 
separations (i.e. comparable to that in RS Oph). 



Recent observations by iPatat et al.l (120071 ) of the Type la SN 2006X show spectral 



features that may be associated with the interaction of the SN ejecta with multiple, slow 
moving shells of pre-existing material, which have been suggested to arise from a recurrent 
nova progenitor system. Furthermore, this interpretation is aided by an anisotropic den- 
sity distribution in the pre-supernova circumstellar environment, compatible with enhanced 
density in the equatorial plane of the progenitor binary system. 

Finally, it is clear that our current spherically-symmetric hydrodynamic models of the 
evolution of shock structures in RS Ophiuchi will need revision to take account of the ge- 
ometry revealed in these observations. Such modeling is now underway (Vaytet et al., in 
preparation). Future work also includes the analysis of a further planned epoch of HST ob- 
servations to explore more details of remant evolution and which may in addition elucidate 
the nature of the more extended features seen in our day 155 observations. In a future paper, 
we hope to combine all our HST observations with contemporaneous high resolution optical 
spectroscopy to enable us to produce the most realistic models to-date of the true geometry 
of the remnant. 
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Object 


Filter 


Exp Time (s) 


RS Oph 


F343N 


240 


RS Oph 


F502N 


180 


RS Oph 


F658N 


40 


HD166215 


F343N 


498 


HD166215 


F502N 


40 


HD166215 


F658N 


16 



Table 1: Summary of total exposure times per filter on both target objects observed with 
the HST ACS/HRC in Program 11004 on 2006 July 17. 
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c) [Ne v] deconvolved 



0.1" 



Fig. 1.— HST ACS/HRC narrow band images of RS Oph (a) PSF-subtracted [O ill] A5007A 
image clearly showing extended emission at sub-arcsecond size scales, particularly in the E- 
W direction (note that the PSF-subtracted [Nev] A3426A image also showed clear evidence 
of extended emission); (b) results of deconvolution of the RS Oph PSF-subtracted image 
with the PSF star, again in [Olll] A5007A, showing a double ring structure; (c) as for (b), 
but for [Nev] A3426A and using a PSF generated by Tiny Tim; (d) deeper, larger area 
view of (b), showing an arc-like feature to the E and a southern "blob" of emission (see text 
for further details). Surface brightness is shown on a logarithmic scale on all images, with 
the vertical bar showing linear representation of the color table from to 50% of the peak 
emission. 



d) [Om] deconvolved 

( 



N 




Fig. 2. — Comparison of the [O ill] A5007 A image with a simple emission model, (a) Side-on 
view of the bipolar nebula, with finite shell thickness. Here, the plane of the central binary 
system lies along the "waist" of the remnant. Surface brightness is simply proportional to 
line of sight path length through the shell; (b) Model inclined by 35 ° to the line of sight, 
consistent with the inclination of the orbit from orbital solutions foobrzycka fc Kenyon 



1994 ); (c) Resulting model with resolution degraded and image rotated 3° anticlockwise 



from an E-W axis to match that of the ACS/HRC [Nev] A3426A image shown in (d) for 
comparison. 



